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I.  SURVEY  OF  ABSORPTION  PHENOMENA  IN  THE  NEAR  INFRARED  REGION 

A.  Attenuation  of  Electromagnetic  Waves  by  a  Dielectric 
Medium:  A  Simplified  Dispersion  Theory 

We  would  like  to  describe,  by  the  use  of  a  simple  model 
the  absorption  of  electromagnetic  radiation  by  non-conducting 
matter.  We  consider  the  material  as  consisting  of  a  collection 
of  N  identical  one-dimensional  oscillators  of  mass  m,  charge  e 
and  elastic  constant  q,  per  unit  volume.  Free  oscillations  are 
described  by 

mvi  +  m<i)02u  =  0  (1) 

in  which  u  is  the  displacement  of  the  particle  from  its  equi¬ 
librium  position,  m  is  the  mass  and  uQ2  =  q/m  is  the  reso¬ 
nant  frequency.  (The  charge  does  not  appear  at  this  point.) 

If  the  oscillator  is  coupled  to  a  loss  mechanism,  a  damping 
term  must  be  added  to  the  above  equation  to  account  for  the 
energy  loss.  The  resulting  damped  oscillator  equation  is 

mu  +  myu  +  mu)2u  =  0  (2) 

o 

in  which  y  is  the  damping  constant. 

The  resonance  frequency  u>0  is  different  from  u)Q .  The  inter¬ 
action  of  the  incident  electromagnetic  field  with  the  oscillator 
is  assumed  to  be  as  follows:  The  periodic  electric  field  of  the 
incident  wave  exerts  a  force  on  the  damped  oscillator  and  drives 
it  at  a  frequency  w.  We  describe  this  motion  with  the  forced 
oscillator  equation 

mvi  +  myu  +  mioQ2u  =  Force  =  eEQe  lwt  (3) 

where  EQ  is  the  maximum  amplitude  of  the  incident  electromagnetic 
wave  oscillating  in  the  u  direction  at  the  site  of  the  oscilla¬ 
tor. 
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We  can  solve  the  above  equation  by  substituting  a  trial 
solution  of  the  form 

u(t)  =  uoe'iut  (4) 

into  equation  (3) .  The  result  is 

eE  /m  .  ^  . 

u  ( t )  =  - - -  e_lwt  =  - e/m -  £(t)  (5) 

to  2  —  oj 2  —  ivui  to  2  —  to 2  —  iY(*> 

0  0 

Thus,  the  amplitude  u  depends  on  the  parameters  of  the  oscil¬ 
lator  (e,m,u>o,Y, )  as  well  as  on  the  strength  of  the  incident 
electromagnetic  field. 

So  far  we  have  considered  just  one  microscopic  process. 

The  description  of  electromagnetic  effects  by  Maxwell's  equations 
uses  macroscopic  quantities  and  material  constants  such  as  the 
dielectric  constant  e  and  permeability  y.  The  connection  be¬ 
tween  the  microscopic  parameters  of  the  oscillators  and  e  and  y 
is  of  interest  here. 

Since  the  amplitude  u,  derived  above,  is  the  displacement 
of  one  oscillator  of  mass  m  and  charge  e,  eu  is  the  induced 
dipole  moment.  The  contribution  of  all  N  induced  dipole  moments 
per  unit  volume  in  the  material  is  the  polarization 

P  =  Neu  (6) 

Introducing  u  from  equation  (5)  above  gives 


m  u)  2  —  u>2  —  iu)Y 

o 


It  should  be  mentioned  that  in  general  E  is  a  function  of  both, 
position  (x)  and  time  (t) . 

We  now  wish  to  solve  the  wave  equation  in  the  medium  mod¬ 
elled  by  the  damped  oscillators.  From  equation  (A-10)  in  appen¬ 
dix  A,  for  a  one-dimensional  dielectric  medium  with  no  conduc¬ 
tivity,  the  appropriate  wave  equation  is 
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3  2E 


3  2P 


(8) 


3  2E 


=  e  y 
2  o  o 


3  x  u  “  3t 

where  we  have  used  the  identity 

1 


—  +  Un 

2  0 


3  t‘ 


=  ye 
o  o 


(9) 


The  second  term  on  the  right  hand  side  of  equation  (8)  is  called 
a  source  term.  It  accounts  for  the  effect  of  polarization 
charges  in  the  medium  as  discussed  in  appendix  A.  We  now  sub¬ 
stitute  for  P  and  E  in  equation  (8)  displaying  explicitly  the 
spatial  and  temporal  dependence  of  the  electric  field,  i.e. 


E  (x , t)  =  EQe 


i  (Kx-ut) 


(10) 


where  K  = 
suit  is 


or 


2tt 


is  the  magnitude  of  the  wave  vector.  The  re- 


-K2  =  -  e  y  u)2  -  y 
0  0  0 


Ne' 


(ui  2  -  w2  -  icoy 


co 


K  ‘ 


Ne2/m£ 


II  + 


U) 


(i)Q2  -  co2  -  icoy, 


(11) 


again  using  equation  (9).  We  see  that  K2  is  complex.  Since  K  is 


Cl) 


related  to  the  index  of  refraction  by  K  = 

r 

the  complex  index  of  refraction  n*=n+ik. 

Then 

K* 2  _  (n  +  ik) 2  _  n2  -  k2  +  2ink 


n,  we  introduce 


(12) 


From  equations  (11)  and  (12)  we  obtain  the  real  and  imaginary 
parts  of  K*2  as 


n2  -  k2  =  1  + 


Ne2  (co  2  -  u)2) 
0 


eQm  [  (u>o 2  -  u)2 )  2  +  ury 2 ) 


(13) 
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r 


and 


2nk  = 


Ne‘ 


e  m 

o 


toy 


(u>  2  -  to2 )  2  +  u)2y‘ 


(14) 


We  call  n  the  refraction  coefficient  and  k  the  extinction  co¬ 
efficient. 

We  have  assumed  above  that  the  driving  electrical  field  at 
the  site  of  the  oscillator  is  equal  to  the  field  of  the  electro¬ 
magnetic  wave  outside  of  the  material.  This  assumption  implies 
that  the  density  of  oscillators  in  the  material  is  low,  somewhat 
as  it  would  be  in  a  gas.  Conversely,  for  a  dense  distribution 
of  oscillators  we  would  have  to  add  the  influence  of  the  sur¬ 
rounding  material,  the  so-called  Lorentz  contribution.  This 
latter  case  will  not  be  considered  here. 

For  such  a  low  density  medium  the  refractive  index  must  be 
close  to  1,  the  value  in  free  space.  As  a  consequence  k  is 
small  compared  to  1  and  we  have  for  n  and  k,  using 
/I  +  x-  -1  +  | 


n  =  1  + 


Ne‘ 


(to  2  -  U)2) 


2e  m 
o 


(to  2  -  to 2 )  2  +  to2y: 


and 


k  = 


Ne" 


2e  m 
o 


J£L 


(u>o  2  -  to2 )  2  +  to 2 y 2 


(15) 


(16) 


It  is  instructive  to  plot  n(oo)  and  k  (to)  and  to  consider  three 

specific  frequencies,  to  =  0 ,  to  =  to  and  to  =  °°.  To  facilitate 

c 

this  discussion,  we  introduce  the  frequency-like  term  u>p=Ne2/eom. 
Then  equations  (15)  and  (16)  become  respectively 


n  (to)  =  1  +  tt  to  2 

2  p 
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2  2 
UJ  -  (A)* 

0 


(toQ2  -  to 2 )  2  + 


2  2 
Y  (A)4 


(17) 


Consider  first  the  refractive  index,  n(co)  plotted  in  Fig¬ 
ure  1 . 


Fig.  1.  Frequency  dependence  of  refractive  index  for  a 
dielectric  described  by  a  simplified  dispersion 
mode 1 . 


For  a)  =  0,  and  y  =  0  (corresponding  to  a  lossless  medium, 
n(0)  =  1  +  h  (a)p/<i>o)  2  =  1  +  (Ne2/2eQq)  .  Numerically,  n(0) 

should  correspond  to  expected  static  values  for  low  density 
media,  i.e.,  n(0)~l.  To  see  this,  we  choose  reasonable  values 

of  N  and  q  and  evaluate  (u>p/ioo)2,  i.e.. 


n  J 


(IQ30  m~ 3 )  (1.6  x  10~19coul)2 
(8.85  x  10-12  farad  m" 1 ) (10 0Nm_ 1 ) 


3  x  10 


U) 


Thus  we  see  that  n(0)  -  1  as  expected.  For  uj  -*•  °°,  n  (<*>)-*- 1 
which  is  the  expected  high  frequency  limit  (for  example,  in  the 
x-ray  region).  For  00  =  (and  y  =  0)  ,  n(w)  is  singular.  This 
situation  is  never  realized  for  real  systems,  however,  since  some 
loss  mechanism  is  always  present.  The  effect  of  the  non-zero 
damping  term  on  n(u))  is  shown  as  a  dashed  line  in  Figure  1.  For 
this  latter  case,  when  to  ~  uQ,  n(to)  >  1  and  is  increasing  as  to 
increases.  This  is  normal  dispersion  behavior  and  is  illustrated 
when  white  light  passes  through  a  glass  prism  with  violet  light 
bending  more  than  red  light.  Near  resonance  (to  =  <oQ)  ,  however, 
n  (to)  decreases  as  to  increases  giving  rise  to  anomalous  dispersion. 
This  effect  would  be  observed  as  the  opposite  as  that  just  des¬ 
cribed,  i.e.  red  light  is  bent  more  than  blue  in  traversing  the 
glass  prism. 

We  now  consider  the  extinction  coefficient  k(<o).  From 
equation  (18),  k(0)  =  0  and  k  (“)-*■  0.  At  resonance  M(oq)  has  a 
maxima,  the  width  depending  upon  y,  the  damping  term.  Clearly 
for  y  =  0,  k(ui)  =  0  for  all  u>.  k(oi)  is  plotted  in  Figure  2. 


Fig.  2.  Frequency  dependence  of  extinction  coefficient  for 
a  dielectric  described  by  a  simplified  dispersion 
model. 
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It  should  be  emphasized  that  the  optical  behavior  just 
described  was  based  upon  a  very  elementary  model  composed  of  N 
identical  oscillators.  Real  materials  require  more  realistic 
models  to  accurately  predict  or  describe  observed  measurements. 
As  an  example  of  a  first  attempt  to  extend  the  above  model,  we 
could  relax  the  restriction  of  identical  oscillators  and  assume 
that  there  are  a  certain  fraction  f ^  of  particles  having  reso¬ 


nant  frequencies 
becomes 


u>  .  =  /  q  .  /m. 
0] 


Then,  for  example,  equation  (11) 


K* 2  = 


OJ 


E 


Ne‘ 


e  m 
o 


f . 


(a>02  -  Ui2  -  iYjOj) 


-] 


(11') 


with  corresponding  expressions  for  n2-k2  and  2nk  following,  as 
before. 

Absorption  Coefficient  a 

We  now  consider  the  effect  of  a  complex  wavenumber  K*  on 
the  propagation  of  a  plane  electromagnetic  wave  in  a  dielectric 
medium.  Equation  (10)  is,  using  K*  =  ^n*=  j^{n+ik)  n+i  ^  k 


i  [£(n+ik)x-u>t] 


E (x, t )  =  EQe 


=  e 


0),  .  ,0)  .  . 

-  =kx  i  (^x-ut) 

u  E  e 

o 


(19) 


Since  the  intensity  I  of  the  radiation  is  the  square  of  the 
electric  field,  we  get 


I  =  I E | 2  =  EE*  =  I E  _ | -  e 


-2(£  k ) x 

2  -  C 


(20) 


Letting  I  s*  lE0l2  and 


a  =  2a  -  2  £  k 


(21) 
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we  get  Beer ' s  Law 

I(x)  =  Ioe  ~ax 


(22) 


Thus  a  tells  us  how  much  a  layer  of  medium  x  units  thick 
attenuates  a  beam  of  radiation  of  intensity  I  incident  upon 
it.  In  particular,  a  ,  the  absorption  coefficient,  measured 
in  units  of  reciprocal  centimeters,  is  such  that  in  a  distance 
(a)"1  the  intensity  falls  from  IQ  to  I (=  IQ/e) .  Experimentally 
a  is  measured  directly.  (See  discussion  below.) 

Using  equation  (21)  and  the  well-known  relations  connecting 
the  various  wave  parameters. 


u>  =  2vf  ,  C  =  Xf 


(23) 


(C  =  speed  of  light  in  vacuo  =  3  x  10 8  m/s;  X  =  wavelength  of 
radiation)  we  see  that 


and 


(24) 


a  is  the  quantity  usually  called  the  absorption  coefficient. 
Equation  (24)  shows  that  a  is  sufficient  to  determine  the  ex¬ 
tinction  coefficient  k.  Note  that  in  entering  a  medium  of  re¬ 
fractive  index  n,  a  plane  wave  (1)  slows  down  by  a  factor  of 
1/n  (since  V  =  C/n) ,  (2)  the  wavelength  decreases  by  a  factor  of 
1/n,  (Xn=  XQ/n;  Xq  =  free  space  wavelength),  but  (3)  the  fre¬ 
quency  remains  constant,  that  is 

“  -  2„f  -  2,  !£-  =  2,  -  a  (25) 


B.  Attenuation  of  Electromagnetic  Waves  by  a  Conducting 
Medium 


Let  us  compare  a  conducting  medium  with  the  dielectric 
medium  previously  considered.  In  a  dielectric  medium,  the 
bound  charges,  represented  as  oscillators,  form  induced  electric 
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dipole  when  exposed  to  incident  electromagnetic  wave.  In  a 
conductor,  however,  this  polarization  mechanism  is  missing  since 
the  electrons  are  not  bound  by  a  restoring  force,  but  are  con¬ 
sidered  to  be  almost  free,  except  at  the  boundaries.  The  inci¬ 
dent  electromagnetic  wave  now  does  not  produce  oscillations  of 
bound  electrons,  but  does  produce  microscopic  currents.  We 
therefore  expect,  in  the  wave  equation,  a  current  density  re¬ 
lated  term.  From  equation  (A-10)  of  appendix  A  we  obtain,  for 
a  one-dimensional  conducting  medium,  a  wave  equation  of  the  form 


=  y  e  J-fS-  +  y  JlL 

3  x2  0  0  3  t2  0  3  t 


(26) 


where  j  is  the  current  density  and  is  related  to  the  velocity 
of  the  electrons  by  j  =  Nev  or  to  the  applied  instantaneous 
electric  field  by  j  =  aE,  a  being  the  electrical  conductivity. 
Again  our  objective  is  to  correlate  the  conductors  microscopic 
parameters  with  the  measured  optical  "constants". 

If  we  try  to  apply  the  damped  oscillator  equation  (with  no 
restoring  force)  we  would  have 

mu  +  myu  =  eEQE  (27) 

The  general  solution  of  this  equation  is  the  sum  of  the  solution 
of  the  homogeneous  and  inhomogeneous  equations.  For  the  homo¬ 
geneous  equation  we  have 

mu  +  myu  =0  (28) 


We  try  as  a  solution 


u  =  e 


t 

T 


and  obtain  y  =  1.  The  damping  constant,  y,  is  inversely  pro¬ 
portional  to  a  decay  time  r  measured  from  an  arbitrary  zero 
point  of  time.  This  time  is  short  (typically  ~  10  1 3  sec)  and 
we  can  neglect  this  contribution. 
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Considering  the  inhomogeneous  equation,  we  observe  that  we 
could  rewrite  it  as 

-iwt 


mv  +  myv  =  eEQe 


(29) 


since  Ci  =  v  and  vi  =  v.  This  is  advantageous  since  we  would 

like  to  correlate  j  =  Nev  to  the  microscopic  process.  A  trial 

.  .  —  iu)t 

solution  v  =  vQe  gives  us 


v  = 


e 

m 


-iu>  +  y 


(30) 


and  for  j  we  obtain 


3  = 


Ne" 


m 


-iu)  +  y 


(31) 


iKx 


Considering  that  j  now  also  depends  on  x  as  e  we  sub¬ 
stitute  for  j  from  equation  (31)  and  E  =  E^'*'  ^Kx  a)t^  into  the 
wave  equation  (26)  ,  and  obtain 


-K2  =  "  W0E0w2  +  UQ 


Ne' 


-10) 


m 


-io>+Y 


or 


k2  =  -ai-  fi  -  Ss 
r2  I 


Let  us  again  set 


K* 2  _  (n  +  ik) 


e  m  io)~y  a) 
o 

2  + 


(32) 


which,  when  substituted 


u)‘  C1 

into  equation  (32)  and  separated  into  real  and  imaginary  parts, 
gives 

Ne2/e  m 

n2  -  Ic2  =  1 - 2 —  (33) 


u)2  +  Y2 


tFor  the  conducting  medium,  the  optical  parameters  will  be 
written  as  n  and  Jc  to  avoid  confusion  with  the  dielectric 
medium  parameters. 
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dliu 


(34) 


Ne2/e  m  y 
2  n  k  =  - 2 —  _ 

u)2  +  y2  to 

These  expressions  are  different  from  the  ones  we  previously 
obtained  for  the  dielectric.  For  this  case,  there  is  no  re¬ 
sonant  frequency  involved  and,  consequently,  no  singularity  pre¬ 
sent  when  y=0»  except  at  co=0.  For  to-*-°°  we  have  k-*0  and  n+1;  which 
means,  for  example  that  x-rays  can  penetrate  metals! 

We  again  use  the  abbreviation  Ne2/eQm  =  to^2  and  plot  n  and 
k  as  functions  of  u).  From  equations  (33)  and  (34)  we  see  that 
for  y^O  and  to-»-0,  n2-k2  is  negative  and  increasing,  that  is 
j  n  |  <  |k|  and  k  has  large  values,  while  for  large  to  we  have  n=l 
and  k  =  0.  From  equation  (33)  we  obtain  n=k  for  00=10^  if  y-0. 

From  the  equation  for  2nk  we  conclude  that  n  is  smaller  than  1. 
Therefore  we  have  qualitatively  for  n  and  k,  the  behavior 
sketched  in  Figure  3. 


The  region  where  k  becomes  small  and  n  approaches  1  is  the 
region  for  which  metals  are  transparent  to  electromagnetic  waves. 
For  large  k,  metals  are  opaque  and  the  metal  reflects  strongly. 
The  wave  can  penetrate  only  a  small  distance  into  the  metal  be¬ 
fore  it  is  completely  attenuated.  An  expression  for  the  pene¬ 
tration  depth  is  obtained  by  the  introduction  of  n  and  k  into  the 
expression  for  the  electromagnetic  wave,  equation  (10),  i.e. 

i[£(n+ik)x-u>t] 

E  =  E  e  u 
o 

-k£x  i [^nx-ut] 

=  EQe  c  e  ^  (35) 

Here  2k  ^  is  the  absorption  coefficient  a.  The  depth  x=6  for 
which  the  wave  is  attenuated  to  e~1  of  its  incident  value  is 
called  the  skin  depth  and  is  obtained  from  ^  k6  =  1  or 

6  =  — —  =  —  (36) 

leu,  a 

The  optical  behavior  of  the  metal,  in  terms  of  the  above  model, 
depends  upon  the  magnitudes  of  id,  and  y  (or  t,  the  relaxa¬ 
tion  time).  For  metals,  typical  values  are  t~10" 1 3 s  (y-10 1 3 s- 1 ) 
and  Up^lO 1 5s~ 1 .  In  the  near  infrared  region  ojsIO^s-1  so  that 
Y=a»<u)p.  Here  k>n  and  for  large  k,  we  see  that  6  will  be  small. 

Note 

Frequently  in  the  literature  there  is  a  different  definition  used 
for  relationships  between  K2  and  n  and  k,  that  is 

(K*)  =  [n  (1  +  ik)]2 

C2 


=  — ^ —  [n2  (1  -  E2)  +  2i  n  1c]  (37) 

C2 
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(38) 


For  the  real  and  imaginary  part  we  would  then  have 


n2 (1  -  k2)  =  1  - 


0) 


u)2  +  y2 


and 


2n  k  = 


oj 2  +  y2  u> 


(39) 


C.  Absorption  of  Electromagnetic  Waves  by  Molecules 

1.  Normal  Vibrations 

In  considering  the  dynamics  of  a  molecule  composed  of 
N  atoms,  3N  coordinates  are  required  to  describe  its  trans¬ 
lational,  rotational  and  vibrational  motion.  If  we  exclude 
linear  and  diatomic  molecules  from  our  consideration,  3 
coordinates  are  used  each  for  the  description  of  the  trans- 
latory  and  rotatory  motion  of  the  entire  molecule.  The  re¬ 
maining  3N-6  coordinates  describe  the  vibrational  motion  of 
the  molecule,  that  is  we  have  3N-6  vibrational  degrees  of 
freedom.  (Linear  and  diatomic  molecules  have  only  3N-5.) 
For  the  determination  of  the  vibrational  frequencies  it  is 
necessary  to  know  the  kinetic  as  well  as  the  potential 
energy  of  the  atom  in  the  molecular  system.  In  the  Born- 
Oppenheimer  approximation,  one  considers  the  electrons, 
which  are  contributing  to  the  bonding  in  the  molecule,  to 
move  much  faster  than  the  nuclei  (with  or  without  attached 
core  electrons).  The  bonding  electrons  are  responsible  for 
the  potential  energy  of  the  force  field  in  which  the  nuclei 
are  undergoing  their  vibrations. 

The  coordinates  used  for  the  description  of  the  vibra¬ 
tional  problem  may  be  initially  of  various  kinds  such  as 
cartesian  coordinates  or  internal  coordinates.  The  kinetic 
and  potential  energy  is  expressed  in  these  coordinates. 
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The  potential  energy  V  is  expressed  as  a  Taylor  series 
about  the  equilibrium  position  of  the  molecule 


V 


3N-6 


'  3  2V  ) 

.9  q^  qj 


o 


q^j 


(40) 


+  (higher  order  terms) 

where  the  ’ s  are  the  mass  weighted  coordinates  of  the 
various  point  masses  (nuclei  with  and  without  core  elec¬ 
trons).  In  the  harmonic  oscillator  approximation,  terms 
higher  than  second  order  are  neglected.  Since  0V/)q^)o 
is,  by  definition  of  the  equilibrium  position,  equal  to 
zero,  and  VQ  is  a  scale  factor  (constant) ,the  potential 
energy  becomes  a  quadratic  expression  in  q^q j ,  i.e. 

3N-6 

2V  =  ^  kii  qi  qi  (A1) 

i » j=l  13  1  3 


where  the  k^'s  are  the  force  constants.  Using  these 
coordinates,  the  kinetic  energy  is  expressed  as 


2T 


3N-6 

I 


i=l 


(42) 


Then,  Lagrange's  formalism  yields  the  equations  of  motion 
3N-6 

q •  +  l  k . .  q.  =  0  j  =  1,  2,  ...,  3N-6  (43) 

3  i=l  13  1 


A  solution  of  this  eigenvalue  problem,  that  is  a  solu¬ 
tion  of  the  system  of  linear  second  order  differential 
equations,  is  obtained  by  introducing  a  trial  solution  into 
the  above  equations : 


qi  =  Ai  cos  +  ^ 


(44) 
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Here  u)=2irf  where  f  is  the  vibrational  frequency  to  be  deter¬ 
mined,  and  if  is  a  phase  factor.  The  systems  of  differential 
equations  yield  a  set  of  simultaneous,  homogeneous  linear 
algebraic  equations  in  the  3N-6  coordinates. 


3N-6 

l  -  <$  •  . u)2 )  A.  =  0  j  =  1  . 3N-6  (45) 

i=l  ^  1 


This  system  of  equations  can  be  solved  only  if  the  secular 
determinant  vanishes,  i.e. 


*il-  «- 


12 


k  -u)‘ 


"1 , 3N-6 
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l3N-6,1 


l3N-6,2  . 


k3N-6,3N-6_a)< 


=  0 


(46) 


Solution  of  the  polynomial  of  order  3N-6  will  yield  at  most 
3N-6  different  u)2,s.  It  may  be  possible  that  some  of  the 
io2,s  are  the  same.  In  such  a  case,  the  associated  vibra¬ 
tional  frequencies  are  said  to  be  degenerate.  For  each  of 
the  3N-6  values  of  io  we  will  obtain  a  set  of  A^'s  from  the 
system  of  algebraic  equations.  These  sets  of  A^'s  are 
called  eigenvectors.  We  can  use  the  3N-6  eigenvectors  to 
compose  a  transformation  matrix  L  of  the  form 


/ 


'll 


Al,3N-6 


1 


(47) 


^  A3N-6,1 


A3N-6 , 3N-6 


I 
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This  matrix  transforms  the  original  coordinate  system  into 
the  normal  coordinate  system. 


(48) 


The  importance  of  this  transformation  can  best  be  seen  if 
we  go  back  to  the  expression  for  the  kinetic  and  potential 
energy.  Both  are  transformed  by  the  above  matrix  from 
their  dependence  on  the  into  dependence  on  Q.  ,  i.e. 


3N-6 

.  ,  L 

3N-6 

2T  = 

l 

i=l 

V - 5 

>1  Qi 

i=l 

3N-6 

2V  = 

l 

j 

kijqiV 

L  ^3N-6 
i=l 

k.Q.  * 
11 


where  the  k^  are  combinations  of  the  k^^ 
motion  appears  in  normal  coordinates  as 


(49) 

The  equation  of 


Qi  +  ki  Qf  *  0  (50) 

that  is,  the  different  equations  are  no  longer  coupled  to 
one  another  and,  therefore,  each  coordinate  is  associated 
solely  with  one  vibrational  mode  of  the  molecule.  In  nor¬ 
mal  modes  all  the  atoms  vibrate  in  phase,  that  is,  they 
all  pass  through  their  equilibrium  positions  at  the  same 
time. 

The  application  of  quantum  mechanics  to  the  system  of 
N  atoms  does  not  change  the  results  drastically.  The 
Schroedinger  equation  is  solved  in  terms  of  normal  coordi- 
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nates.  The  inodes  of  the  molecular  motion  are  associated 
with  wave  functions  and  eigenvalues.  The  classical  vibra¬ 
tional  frequency  is  the  same  as  the  one  obtained  from 
quantum  mechanics.  However,  selection  rules  are  introduced 
and,  consequently,  the  absorption  intensity  is  modified. 


2.  Selection  Rules  for  a  Vibrating  Molecule 

A  vibrating  molecule  may  absorb  a  light  quantum  hf, 
where  f  corresponds  to  the  energy  difference  between  the 
initial  and  exited  (finite)  state.  The  Bohr  frequency  con¬ 


dition  is 


E,  -  E. 
f  1 


where  h  =  Planck's  constant.  The  selection  rules  are  ob¬ 
tained  by  expanding  the  electric  dipole  moment  y  of  the 
molecule  in  terms  of  normal  coordinates,  i.e. 


3N-6 


y  =  u0  + 


N-0  r  1 

7  — —  Qv  +  (higher  order  terms)  (52) 

K=1  *3  Q  J  o  K 


Again,  in  the  harmonic  oscillator  approximation,  the  higher 
order  terms  are  neglected  for  the  dipole  moment  matrix 
element  /  lj^dx-  We  then  have 


3N- 

U±Ujfdr  +  l 


6MA 

1  l  9  qkJ  J  1 


VfdT  ■ 


i|>.  may  be  written  as  a  product  of  the  initial  states  of  all 


vibrations  3N-6 

=  n  i|» 


and,  if  only  the  j  state  is 


singly  exited,  may  be  written  as 

3N-6 

^  .  n  ip  „ 
f  K^j  0K 
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Introduction  of  the  wave  function  into  the  expression  for  the 
matrix  element  of  the  dipole  moment  yields 


The  first  integral  vanishes  since  y  is  a  constant  and 

dx  =  0  from  the  orthogonality  condition.  In  the 

second  integral  all  terms  in  the  sum  of  the  type  /i|i  „  Q„  ip  v 

OK  K  OK 

are  zero  by  symmetry.  The  only  non- zero  term  is 


dx 


(56) 


would  vanish  or  if  we  would 
only  one  vibration  to  be 


This  term  would  be  zero  if 
not  have  initially  chosen 
singly  excited. 

The  general  selection  rule  for  harmonic  oscillator  vibra¬ 
tional  transition  is 


Aj  =  ±  1  (57) 

The  permanent  dipole  moment  has  no  influence  on  the  vibrational 
transition,  but  the  change  in  the  dipole  moment  with  respect  to 
the  vibrational  coordinate  (QK)  is  important  and  must  be  non¬ 
zero.  The  harmonic  oscillator  approximation  is  sufficient  to 
explain  the  most  intense  transitions.  By  considering  the  in¬ 
fluence  of  the  higher  order  terms  in  equations  (40)  and  (55)  on 
the  dipole  moment  y,  weaker,  previously-forbidden,  transitions 
may  be  predicted. 
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3 .  The  Concept  of  Group  Frequencies 

Reviewing  the  near  infrared  absorption  spectra  of  a  large 
number  of  chemical  compounds,  one  observes  that  absorption 
bands  at  certain  frequencies  are  related  to  the  occurence  of 
certain  molecular  groups  in  the  molecule.  For  example,  the 
bands  around  the  3000  cm-1  region  are  related  to  C-H  stretch¬ 
ing  vibrations.  A  careful  study  of  this  phenomenon  has  re¬ 
vealed  that  the  relationship  between  observed  bands  and  molec¬ 
ular  groups  is  almost  unique  and  therefore  one  speaks  of  the 
fingerprint  region  if  the  analytical  aspect  of  the  near  infra¬ 
red  region  is  considered.  The  correspondence  between  normal 
vibrations  and  group  frequencies  is,  for  most  of  the  cases, 
simple.  Let  us  consider  the  normal  mode  of  a  certain  fre¬ 
quency  which  is  also  listed  as  group  frequency.  The  atoms  be¬ 
longing  to  the  group  frequency  vibrations  are  undergoing  large 
amplitude  vibrations  in  the  normal  mode  compared  to  all  the 
other  atoms.  From  the  point  of  view  of  interaction  with  the 
electromagnetic  field  the  effect  of  all  the  other  atoms  may  be 
neglected. 

D.  Lattice  Vibrations  in  Dielectric  Media 
1 .  Bulk  Material 

The  periodic  arrangement  of  molecules  in  a  crystal  lat¬ 
tice  leads  to  new  types  of  absorption  phenomena  which  are 
observed  in  addition  to  the  absorption  due  to  the  internal 
vibrations.  For  most  molecules  in  crystal  form,  lattice  vi¬ 
bration  absorption  is  observed  in  the  middle  and  far-infrared 
regions.  Only  for  small  molecules  may  we  observe  this  type 
of  absorption  in  the  near  infrared.  Consequently,  the  one 
dimensional  linear  chain  may  well  serve  as  a  model  for  the 
discussion  of  this  phenomena.  In  this  model,  the  unit  cell 
contains  two  atoms  with  different  effective  charges  so  that 
a  contracting  force  exists  between  an  atom  and  its  nearest 


neighbors.  One  can  see  immediately  the  occurrence  of  two 
fundamental  modes ,  one  in  which  the  two  atoms  of  the  unit 
cell  move  in  opposite  directions  and  one  in  which  they 
move  in  the  same  direction.  The  former  is  called  the  longi¬ 
tudinal  optical  mode,  with  frequency  <i>T  ,  whereas  the  latter 
is  the  longitudinal  acoustical  mode,  with  zero  frequency.  A 
more  detailed  analysis  leads  to  the  dispersion  relation  of 
the  modes,  that  is,  the  relationship  between  the  frequency 
and  wave  vector  for  all  possible  modes  in  the  Brillouin  zone. 

If  we  extend  these  considerations  to  three  dimensions,  we 
find  that  there  are  two  transvers  optical  and  two  transverse 
acoustical  mode  branches  possible,  in  addition  to  the  longi¬ 
tudinal  one.  Of  all  these  modes,  we  are  only  interested  in 
those  which  lead  to  strong  optical  absorption.  For  inter¬ 
actions  with  an  electromagnetic  field  to  occur,  we  need,  as 
usual,  a  component  of  the  electric  field  in  the  direction  of 
the  change  of  the  dipole  moment  of  the  chain.  Since  the 
transverse  modes  are  degenerate ,  we  have  only  two  branches 
of  modes  to  consider.  Selection  rules  predict  that  fre¬ 
quencies  belonging  to  zero  wave  vectors  are  the  fundamental 
frequencies.  Difference  and  summation  frequencies  may  occur 
at  other  non-zero  wave  vectors  with  far  lesser  intensity. 

The  fundamental  frequencies  can  easily  be  visualized  by  con¬ 
sidering  the  chain  as  composed  of  two  sublattices  of  equal 
atoms  having  sublattices  oscillating  rigidly  against  one 
another.  This  results  in  a  large  change  in  the  dipole  mo¬ 
ment  and,  consequently,  in  a  large  absorption.  The  absorp¬ 
tion  corresponding  to  this  fundamental  vibration  is  called, 
in  ionic  crystals,  the  reststrahlen  band.  The  names  re¬ 
late  to  the  fact  that  in  very  strong  absorption,  the  incident 
wave  is  attenuated  in  such  a  small  surface  layer  that  the 
surface  reflects  most  of  the  radiation,  as  is  found  for  a 
metal. 
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Theoretically  derived  quantitative  expressions  for  the 
absorption  coefficient  are  difficult  to  obtain  since  it  is 
necessary  to  know  the  wave  functions  of  the  ground  and 
excited  state  in  order  to  calculate  the  change  in  the  di¬ 
pole  moment.  Experimentally,  one  proceeds  as  discussed  at 
the  beginning  of  this  section,  and  determines  k,  from 
which  a  can  be  obtained. 


2.  Dielectric  Powders 

We  have  just  seen  that  the  fundamental  absorption  pro¬ 
cess  in  ionic  type  crystals  involves  the  transverse  optical 
mode.  In  special  cases,  it  may  also  involve  the  longitudi¬ 
nal  optical  mode.  If  the  diameter  of  the  crystal  becomes 
very  small ,  the  frequency  of  the  fundamental  vibration 
shifts  to  higher  values.  It  has  been  observed,  e.g.  in 
MgO  micro-crystals  of  radius  50A  to  800A,  that  the  lattice 
absorption  peak  shifts  up  by  10  cm  The  reason  for  this 
shift  can  be  seen  by  considering  the  number  of  atoms  on  the 
surface  of  the  crystals  (having  neighbors  partially  surround¬ 
ing  them)  compared  to  the  number  of  atoms  in  the  interior  of 
the  crystal  (having  neighbors  completely  surrounding  them) . 

For  a  particle  radius  of  50A,  about  25%  of  the  ions  are  on 
the  surface,  whereas  for  a  radius  of  380A  there  are  only  3% 
on  the  surface.  Since  the  shift  is  to  higher  frequencies 
with  decreasing  size,  one  may  conclude  that  the  lattice 
spacing  on  the  surface  is  smaller.  This  may  indeed  be  proven. 
As  a  consequence  of  this  finding,  one  expects  a  strong  in¬ 
fluence  on  the  lattice  vibrational  frequency  of  the  surround¬ 
ing  matrix  due  to  the  polarization  charges.  One  finds,  how¬ 
ever,  a  different  frequency  oiF  which  is  related  to  the  trans¬ 
verse  optical  mode  frequency  by 
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where  eQ  and  e#  are  the  low  and  high  frequency  dielectric 

constants  and  e  the  dielectric  constant  of  the  surrounding 
in 

medium.  The  frequency  shift  due  to  this  effect  can  be  re¬ 
markable.  In  BaTiO.  ,  for  example,  the  lattice  vibrational 

’  -1 

frequency  of  the  bulk  material  is  at  33  cm  ,  whereas  for 

small  spheres  in  powder  form  in  a  surrounding  medium  with 

e  =  1 ,  one  observes  a  value  of  400  cm-1, 
m 

Experiments  on  small  particles  are  hampered  by  the  scat¬ 
tering  effects  of  the  particles.  Quantitative  experimental 
values  can  be  obtained  via  the  dielectric  constant  and  from 
that,  the  determination  of  a  may  be  made. 

E.  Absorption  of  Electromagnetic  Waves  by  Semiconductors 

For  absorption  by  semiconductors  in  the  near  infrared 
spectral  region  we  consider:  (a)  Transitions  of  electrons 
across  the  band  gap,  (b)  free  carrier  absorption  and  inter¬ 
band  absorption,  and  (c)  multiphonon  absorption. 

1.  Band  Gap  Absorption 

The  strongest  absorption  in  the  near  infrared  is  due  to 
the  transition  of  electrons  from  the  valence  band  to  the  con¬ 
duction  band.  These  transitions  are  only  possible  if  the 
energy  of  the  absorbed  photon  is  comparable  or  larger  than 
the  band  gap  energy,  that  is  hf  i  E  .  For  photons  with  hf 
smaller  them  Eg  these  transitions  are  not  possible.  A  plot 
of  the  absorption  coefficient  a(f)  against  the  photon  energy 
shows  therefore  a  strong  increase  of  absorption  near  hf=Eg; 
this  is  called  the  band  edge  absorption.  For  hf  <  Eg  the 
material  is  relatively  transparent.  The  absorption  in  the 
hf  >  Eg  region  is  relatively  strong;  a  has  values  of  the 
order  of  10 ‘‘/cm. 

The  transition  from  the  valence  to  the  conduction  band 
can  be  a  "direct"  or  an"indirect"transition,  depending  on 
the  structure  of  the  bands.  The  energy  difference  Eg 


30 


between  the  bands  depends  in  general  on  the  crystal  struc¬ 
ture.  Direct  transitions  occur  if  the  minimum  of  the  con¬ 
duction  band  and  the  maximum  of  the  valence  band  both  are 
situated  at  the  same  K  value  (e.g.  at  K  *  0).  The  electron 
then  makes  a  "vertical  transition"  without  change  in  mo¬ 
mentum  (K  is  conserved)  and  the  energy  difference  is  E  =  hf. 

9 

An  indirect  transition  occurs  if  the  minimum  of  the  conduc¬ 


tion  band  and  the  maximum  of  the  valence  band  are  not  situ¬ 


ated  at  the  same  point  in  K  space.  In  order  to  conserve  mo¬ 
mentum  a  phonon  is  either  absorbed  or  emitted.  Since  the 
photon  momentum  is  small  in  the  energy  region  under  consid¬ 
eration,  the  phonon  makes  up  for  the  change  in  K.  On  the 
other  hand,  the  phonon  energy  is  small  compared  to  the  energy 
of  the  absorbed  photon  so  that  hf  =  E^.  The  absorption  pro¬ 
bability  for  an  indirect  process  is  smaller  that  that  for  a 


direct  process. 


Let  us  assume  a  parabolic  relationship  between  E  and 
K  at  the  maximum  of  the  valence  band  and  at  the  minimum 
of  the  conduction  band  where  the  direct  transition  occurs. 
The  absorption  coefficient  is  then  found  to  be: 


q  is  the  charge  of  the  electron,  m^*  and  mg*  are  the 
effective  masses  of  holes  and  electrons  respectively,  n 
is  the  real  refractive  index,  C  is  the  speed  of  light, 
and  h  is  Planck's  constant.  For  n  =  4  and  m^*  =  me*  =  mg  , 
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(61) 


where  me  is  the  free  electron  rest  mass,  we  have 
a (f )  =  2  x  10“  (hf  -  E  )* 

g 

a(f)  is  given  in  cm-1  if  h  and  E  are  inserted  in  units  of 

g 

electron  volts. 

(b)  Absorption  coefficient  for  forbidden  direct 
transitions 

Sometimes  direct  transitions  are  forbidden  by  selec¬ 
tion  rules  for  K  =  0  (symmetry) ,  but  allowed  for  K  ?  0. 

For  this  case  the  absorption  coefficient  is  found  to  be: 

a(f)  =  A'  (hf  -  E  )  3/2  (62) 

g 

where 


A'  =  1 

A  3 


q  V  +  me 

n  C  h2  me*  mh*(hf) 


and  for  the  case  discussed  above  [paragraph  (a) ]  we 


have 


ct(f)  =  1.3  x  10‘ 


(hf  -  E  )3/z 


,  cm 


again,  hf  and  Eg  are  in  ev's. 

(c)  Absorption  coefficient  for  indirect  transitions 

The  indirect  transitions  are  two-step  processes  and, 
consequently,  the  interaction  of  photons  with  electrons 
and  phonons  have  to  be  considered.  The  absorption  co¬ 
efficient  has  been  shown  to  be 
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where  rs  the  energy  of  the  phonon.  The  first  term  is 
due  to  a  process  with  phonon  absorption,  the  second  with 
phonon  emission.  The  constant  A  has  a  complicated  ap¬ 
pearance.  In  many  cases  it  is  determined  by  an  experi¬ 
mental  fit  to  the  data. 

2.  Free  Carrier  Absorption 

Free  carrier  absorption  is  of  importance  in  semicon¬ 
ductors  as  it  is  in  conductors.  Absorption  due  to  free  car¬ 
riers  may  be  observed  in  the  near  and  far  infrared.  The  ab¬ 
sorption  of  photons  by  the  electrons  in  the  periodic  field 
of  the  crystal  is  strictly  not  allowed  as  a  consequence  of 
conservation  of  momentum.  However,  a  perturbation  involving 
lattice  vibrations  (phonons)  relaxes  the  restrictions  and  ab 
sorption  of  photons  is  possible.  The  absorption  coefficient 
may  be  calculated  on  the  basis  of  classical  theory  by  assum¬ 
ing  an  oscillating  electron  without  a  restoring  force  and 
using  an  effective  mass  me*  in  place  of  the  electron  rest 

mass  m  .  One  obtains: 
e 


n2  -  k2  -  e  = 


2nku)  = 


Nq2/me*£0 
to2  +  y2 
Nq2/me*e0 
a)2  +  Y2 


where 


and  t  is  the  collision  time.  For  wave¬ 


lengths  up  to  a  few  hundred  micrometers,  the  absorption  co¬ 
efficient  may  be  expressed  as: 


4ir2C3n  e. 


m  2u 


where  X  is  the  wavelength;  q  is  the  electron  charge;  me  is 
the  electron  mass;  ue  is  the  electron  mobility;  c  is  the 
speed  of  light;  n  is  the  refractive  index  and  eQ  is  the 
static  dielectric  constant.  The  order  of  magnitude  of  the 
absorption  coefficient  a(f)  is  10  -  100  cm-1  for  free  car¬ 
rier  absorption  in,  for  example,  n-type  InAs  in  the  3-15ym 
region.  This  value  is  two  orders  of  magnitude  smaller  than 
that  for  band  edge  absorption.  Interband  absorption  and  ab¬ 
sorption  due  to  holes  results  in  absorption  coefficients  of 
the  same  magnitude. 

3.  Multiple  Phonon  Absorption 

The  absorption  due  to  lattice  vibrations  occurs,  for 
most  material,  in  the  middle  and  far  infrared.  For  ionic 
crystals  these  lattice  vibrations  show  up  as  reststrahlen 
bands.  These  bands  are  due  to  fundamental  modes  of  the  lat¬ 
tice  with  K  =  0,  K  being  the  crystal  momentum  vector.  Mul¬ 
tiple  phonon  absorption  of  the  summation  type  can  be  ob¬ 
served  in  the  near  infrared  region.  Absorption  due  to  such 
processes  are  relatively  weak,  because  these  processes  are 
second  order  effects.  They  appear  as  modulation  on  top  of 
the  free  carrier  absorption. 

F.  Absorption  of  Electromagnetic  Waves  by  Metals 

1.  Bulk  Materials 

The  dependence  of  the  optical  constants  on  the  micro¬ 
scopic  parameters  of  bulk  conducting  materials  has  been  dis¬ 
cussed  in  the  first  chapter  on  "simplified  dispersion  theory". 
In  this  simple  approach  the  plasma  frequency  and  relaxation 
time  are  the  characteristic  parameters  for  the  absorption 
process.  This  process  is  similar  to  the  absorption  in  semi¬ 
conductors  as  we  have  just  seen  in  section  I.  E.2  "free 
carrier  absorption" . 
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2.  Metal  Powders 


By  considering  metal  powder  instead  of  the  bulk  mat¬ 
erial  one  anticipates  changes  analogous  to  those  observed 
for  dielectric  bulk  material  and  their  powders.  The  impor¬ 
tant  parameter  is  the  plasma  frequency  and  for  particles  of 
about  10A  diameter,  a  broadening  of  the  plasma  resonance 
peak  is  observed.  Theoretically,  this  broadening  is  des¬ 
cribed  by  a  decreased  electron  mean  free  path  due  to  scat¬ 
tering  at  the  particle  surface.  This  macroscopic  explana¬ 
tion  has  been  more  successful  than  considerations  from  a 
quantum  mechanical  point  of  view.  In  thin  conducting  films 
on  the  other  hand,  one  has  observed  effects  which  are  inter 
preted  as  the  breakdown  of  the  continuous  conduction  band 
into  discreate  states  by  reducing  the  thickness  of  the  film 
from  bulk  material  dimensions.  The  evaluation  of  the  ab¬ 
sorption  coefficient  from  experimental  measurements  pro¬ 
ceeds  as  outlined  before. 

G.  Absorption  of  Electromagnetic  Waves  by  One-Dimensional 
Conductors 

Recently  a  number  of  linear  conducting  materials  have 
been  developed.  These  polymeric  materials  have  metallic  be¬ 
havior  in  the  direction  of  the  chain  and  almost  insulator- like 
properties  in  the  direction  perpendicular  to  the  chain.  An  ex¬ 
ample  of  a  polymeric  one-dimensional  conductor  (not  recognized 
as  such  initially)  is  the  H-sheet  polarizer  invented  by  E.  Land 
in  1938.  This  polarizer  is  made  of  polyvinyl  alcohol  chains 
which  are  aligned  by  uniaxial  stretching  of  the  material.  The 
sheet  is  impregnated  with  iodine  by  immersing  it  into  a  special 
ink.  The  iodine  atoms  or  molecular  clusters  are  believed  to  do¬ 
nate  electrons.  These  electrons  constitute  the  free  charge  car¬ 
riers  involved  in  the  conduction  process  as  in  a  metallic  con¬ 
ductor.  Incident  light  with  its  electric  field  vector  oriented 
parallel  to  the  polymer  chain  is  absorbed  so  that  the  sheet  acts 
as  a  polarizer. 
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As  a  result  of  the  considerable  interest  in  linear  conductors , 
a  number  of  one-dimensional  linear  polymeric  conductors  have  been 
extensively  studied,  e.  g.  (SN)x  and  polyacetylene-iodide.  These 
materials  are  black  and  light-weight,  and  for  polyacetylene-iodide, 
the  absorption  in  the  near  infrared  region  is  broad  and  strong, 
but  without  distinctive  absorption  features.  A  quantitative  ab¬ 
sorption  study  is  underway  in  our  laboratory. 

II.  QUANTITATIVE  ANALYSIS  OF  ABSORBING  MATERIALS 

An  electromagnetic  wave  incident  upon  a  material  will  inter¬ 
act  with  it,  ar.d  have  its  energy  changed,  by  three  types  of  mech¬ 
anisms:  (1)  Reflection,  in  which  a  fraction  of  the  incident  radi¬ 

ation  is  deflected  backwards;  (2)  scattering,  in  which  a  fraction 
of  the  incident  radiation  is  scattered  in  all  directions  (charac¬ 
teristics  depend  upon  relationship  between  particle  size  and  wave¬ 
length  of  radiation,  optical  inhomogeneities,  etc.);  and  (3)  ab¬ 
sorption,  the  general  term  describing  the  conversion  of  radiant 
energy  into  thermal  energy.  Generally,  all  these  processes  are 
present  to  varying  degrees  and,  are  sensitive  to  the  wavelength 
of  the  incident  radiation.  We  are  interested  here  in  the  ab¬ 
sorption  process  and  assume  that  the  beam  incident  normally  upon, 
and  entering  the  medium,  has  a  radiant  power  IQ.  Upon  passing 
through  a  uniform  layer  b  units  thick,  the  radiant  power  is 
attenuated  to  a  value  I (b)  given  by  the  absorption  (Lambert- 
Bouguer)  law 

1(b)  =  I  e  "a<v)b  (69) 

o 

The  product  a(v)b  is  the  optical  density  or  absorbance  and  a(v) 
is  the  (frequency-dependent;  v  =  i-)  absorption  coefficient.  As 
rantioned  in  paragraph  I.  A,  a(v)  is  related  to  the  imaginary 
part  of  the  complex  index  of  refraction  k(v)  =  Im[n(v)*]  by 

a(v)  =  — iy  k(v)  =  4irvk(v)  (70) 
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The  extinction  coefficient  k(v)  is  a  function  of  frequency  [as 
is  n(v)].  Equations  (69)  and  (70)  show  that,  in  principle,  one 
has  only  to  measure  the  total  sample  thickness  b  and  to 

extract  values  of  a(v)  and,  if  desired  k(v),  for  material^  in 
the  form  of  uniform  slabs.  If  the  material  of  interest  is  dis¬ 
solved  in  a  suitable  (i.e. ,  non-absorbing,  non-interacting,  etc.) 
solvent,  Eqn.  (69)is  written  as 


1(b) 


r  “  e  (v)  cb 
I  e 


(71) 


usually  referred  to  as  Beer's  Law.  Here  b  is  the  optical  path 
length  through  the  solution  of  concentration  c  moles/liter 
characterized  by  a  molar  extinction  coefficient  e(v)  measured 
in  units  of  liter/mole-cm.  Clearly  for  a  liquid  film  confined  to 
a  cell  of  thickness  b  cm,  comparison  of  equations  (69), (70)  and 
(71)  gives 

e(v)  =  — vk(v)  (72) 

The  concentration  may  be  obtained  from 


c  = 


(73) 


where  p  is  the  solute  density  (in  gm/mi.)  ,  M  is  the  molecular 
mass  (in  gm/mole)  ,  and  4>  is  volume  fraction  of  solute  [i.e.,  mJ. 
(solute) /liter  (solution)].  Clearly,  in  assuming  Beer's  Law 
to  hold  for  a  solution,  we  must  determine  experimentally  the 
absorbance  A  and  either  e(v)  or  c.  In  routine  analysis  of  mat¬ 
erials,  e(v)  is  usually  known  and  the  concentration  c  is  the 
desired  quantity. 

A  quantity  of  considerable  interest,  especially  as  related 
to  the  study  of  molecular  dynamics,  is  the  integrated  absorption 
intensity,  defined  here  as 


V  2 

a (v)dv 


Vi 


(74) 
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Here  and  are  frequencies  (wave  numbers)  sufficiently  far 

removed  from  the  band  center  that  a(v  )  =  a(v  )  =  0. 

1  2 

Theories  relating  the  microscopic  processes  to  a(v)  are 
mathematically  quite  involved.  The  resonance  absorption  of  a 
molecule  in  the  gaseous  state,  for  example,  depends  on  both  the 
vibrational  frequency,  and  the  frequency  of  the  overall  rotation 
of  the  molecule.  Interactions  between  the  vibrational  and  rota¬ 
tional  motion  occur.  The  analysis,  theoretical  or  empirical, 
should  reveal  the  observed  frequencies.  The  intensities  may  then 
be  calculated  by  assuming,  for  example,  a  Boltzmann  equilibrium 
distribution  of  the  oscillators.  The  effect  of  the  collisions 
between  the  molecules,  that  is,  line  broadening,  might  be  con¬ 
sidered  as  well  as  other  factors  affecting  the  microscopic  pro¬ 
cesses,  and  therefore  the  absorption  coefficient.  Calculations 
are  only  possible  for  the  simplest  of  molecules.  For  larger 
molecules  one  needs  interaction  constants  of  the  vibration- 
rotation  interaction  or  the  vibrational  wave  functions  of  the 
ground  and  first  excited  state.  These  are  difficult  to  calculate 
and  empirically  not  easy  to  determine.  As  an  "over-simplified" 
rule,  however,  in  most  cases  the  absorption  coefficient  is  pro¬ 
portional  to  the  resonance  frequency  and  the  square  of  the  matrix 
element  of  the  dipole  moment.  In  liquids,  the  overall  rotation 
of  the  molecule  is  suppressed  and  the  collisions  between  molecules 
play  a  more  dominant  role.  The  absorption  process  might  be  con¬ 
sidered  as  consisting  of  resonance  absorption  and  orientational 
absorption.  The  former  is  an  absorption  centered  about  the  vi¬ 
brational  resonance  frequency  where  the  rotational  states  do  not 
exist  as  individual  levels,  but  contribute  to  a  broad  absorption 
band.  The  latter  depends  on  the  orientation  (and  reorientation) 
of  permanent  dipoles  in  an  applied  electromagnetic  field  (Debye 
absorption) .  An  empirical  description  of  these  two  processes 
is  of  the  form 

a(t)  =  a  e  t//x  cos  (wt  +  4>)  +  a  e  t'^T 
1  2 
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If  the  constants  a  ,  t,  a  and  <J>  are  determined,  a(v)  may  be  ob- 

12 

tained  through  a  Fourier  transformation.  Here  t  is  a  character¬ 
istic  relaxation  time  associated  with  the  orientation.  Theoret¬ 
ical  calculations  of  a(v)  using  procedures  from  statistical  mech¬ 
anics,  as  well  as  quantum  mechanics,  are  difficult  and  lengthy. 
However,  this  is  a  very  active  field  since  one  can  deduce  the 
correlation  function  and  obtain  an  estimate  of  the  molecular  be¬ 
havior  subsequent  to  an  assumed  initial  configuration.  The  time 
interval  considered  is  of  the  order  of  10" 11 -10“ 12  sec. 

Experimental  Methods 

The  use  of  the  Lambert-Bouguer-Beer  Law  appears  deceptively 
straightforward.  While  the  method  is  not  sophisticated,  there 
are  certain  precautions  that  must  be  considered  so  that  the  data 
yield  true  material  parameters  rather  than  experimental  artifacts. 
We  will  consider  in  limited  detail,  the  more  significant  of  these 
features . 

A.  Liquids 

1.  Deviations  from  Beer's  Law  and  Photometric  Error 

The  attenuation  of  a  beam  of  light  as  it  passes  through 
a  specimen  depends  upon  the  samples  thickness  b,  its  ab¬ 
sorption  coefficient  a  and  the  particular  wavelength  A  or 
frequency  (actually  wavenumber  v  *  A-1).  Ultimately,  the 
attenuation  depends  upon  the  number  of  absorbing  particles 
(molecules)  "seen"  by  the  optical  beam.  Provided  that  the 
number  of  particles  and  the  volume  containing  them  remain 
constant,  the  absorption  coefficient  a(v)  suffices  to  char¬ 
acterize  the  substance.  Such  is  the  case  for  solids  and 
pure  liquids.  For  solutions  or  mixtures  of  liquids,  however, 
this  condition  is  not  fulfilled  and  it  becomes  necessary  to 
relate  the  absorption  to  the  sample  concentration.  (The 
analogous  situation  for  gases  would  require  that  the  effect 
of  pressure  be  considered.)  This  is  done  in  the  case  of 
Beer's  Law,  since  the  absorbance  is  given  by 
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A  =  c  (v)  c  b  (75) 

For  a  substance  to  "obey"  Beer's  Law  then,  the  absorbance 
must  be  proportional  to  the  concentration.  Unfortunately, 
this  is  not  always  found  and  deviations,  both  real,  as  in 
the  case  of  refractive  index  changes  due  to  high  solute  con¬ 
centrations,  and  apparent,  due  to  instrumental  limitations 
and  non-symmetrical  chemical  equilibrium,  are  encountered. 

The  real  deviations  arise  when  the  measured  absorption  is  no 
longer  the  sum  of  the  absorptions  for  all  separate  particles. 
More  specifically  the  relation 

N  N 

A  =2  Ai  -  b  2  Ei<V)Ci  (76) 

i=l  i=l 

where  i  labels  the  species  and  N  is  the  total  number  of 
species  present,  is  not  obeyed.  This  non-ideal  behavior  is 
a  result  of  formation  of  molecular  aggregates  or  clusters  by 
solute-solute  and  solute-solvent  interactions.  This  type  of 
deviation  may  be  used  to  study  molecular  interactions,  and  is 
therefore,  not  without  some  value. 

The  experimental  deviations  arise  because  Beer's  Law  is 
rigorously  correct  only  for  monochromatic  radiation.  The 
finite  slit  width  of  the  instrument  monochromator  allows  a 
band  of  frequencies  to  be  detected  so  that  the  magnitude  of 
the  deviation  depends  on  the  ratio  of  the  slit  width  to  the 
absorption  band  width  observed.  This  effect  has  no  compen¬ 
sating  benefits,  but  must  be  tolerated.  However,  empirical 
determination  of  factors  affecting  the  absorbance  may  still 
be  evaluated  provided  the  aforementioned  instrumental  fac¬ 
tors  are  maintained  constant. 

The  above  considerations  of  deviations  leads  to  some 
practical  rules-of-thumb  for  choosing  the  experimental  para¬ 
meters.  To  eliminate  clustering  effects  in  solutions,  con¬ 
centrations  below  about  .01M  should  be  used.  However,  the 
instruments  ability  to  detect  small  changes  in  radiant  power 
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limits  the  dilution  to  about  10" 7M.  Concommitant  with  the 
desire  for  low  concentrations  is  the  instrument's  need  for 
sufficient  energy  at  the  detector  so  that  photometric  error 
may  be  minimized.  Here  again  a  suitable  balance  must  be 
sought  since  for  a  low  solute  concentration,  a  weak  ab¬ 
sorption  band  results  and  the  instrument  "sees"  a  small 
difference  between  incident  and  transmitted  power,  and  thus 
a  relatively  large  error  in  absorbance.  Alternatively,  for 
a  more  concentrated  solution,  the  absorption  band  is  strong 
and  insufficient  power  reaches  the  detector,  also  giving  in¬ 
accurate  results.  An  appropriate  intermediate  range  of  ab¬ 
sorbance  can  be  obtained  starting  with  Eqn.  (75)  rewritten 
as 

E(u)  -  -is  -  -h  1o?  -I - h  109  T  <77) 

Here  we  have  used  the  identity  A  =  -log  T,  where  T  is  trans¬ 
mission,  usually  associated  with  I/IQ.  Our  aim  is  to  mini¬ 
mize  the  error  in  the  molar  extinction  coefficient  e (v)  by 
minimizing  the  error  in  T.  The  fractional  error  in  e(v)  is 

ds (v)  _  _  log  e  dT  ,_Q. 

e (v)  e (v)bc  T 

where  log  means  log  ("Jin"  is  reserved  for  log  )  log  e  = 

10  ® 

1/2.30  =  .434,  and  dT  is  the  uncertainty  in  T.  To  minimize 
de/e,  differentiate  Eqn.  (78)  and  equate  the  result  to  zero. 
The  result  gives  log  T  +  log  e  =  0,  log  T  =  -.434,  or  TQpt  = 
.368  for  the  optimum  value  of  transmission  to  minimize  in¬ 
strument  photometric  error. 

To  summarize,  optimum  accuracy  in  the  determination  of 
e(v)  is  obtained  when  the  solute  concentration  is  such  that 
(1)  molecular  clustering  is  avoided  and  (2)  the  transmission 
is  in  the  range  of  30-40%  for  absorption  band  minima. 


2.  Cell  Losses 

In  stating  the  Lambert-Bouguer- Beer  Law  [Eqns.  (69)  and 
(71) ]  the  quantity  I  is  taken  as  the  incident  beam  power  and 
I  the  transmitted  power.  Since  a  beam  of  light  is  partially 
reflected  at  the  interfaces  between  two  dielectric  media, 
the  incident  power  arriving  at  the  first  surface  of  a  liquid 
cell  is  not  totally  transmitted  to  the  second  surface,  where 
additional  reflection  also  takes  place.  In  fact,  multiple 
reflections  will  take  place  at  all  four  interfaces  of  a  liquid 
cell  consisting  of  two  (assumed)  flat  windows  separated  by 
an  (assumed)  uniform  space.  Ignoring  the  convergence  of  the 
optical  beam  in  the  sample  compartment  of  most  infrared 
spectrometers,  the  fraction  of  incident  radiation  reflected 
at  normal  incidence  can  be  calculated  from 
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For  an  NaCl  window  (n2  =  1.54)  in  air  (nx  =  1.00), 

Irefl/^Io  =  or  4.5%  of  the  incident  power  is  reflected 

from  the  first  surface.  Clearly,  for  accurate  quantitative 
work,  the  various  cell  losses  must  be  accounted  for.  One 
method  suggested  by  Draegert,  et  al  [J.O.S.A.  5£,  64,(1966)] 
and  used  in  this  reported  work  is  to  consider  the  actual  in¬ 
cident  beam  to  be  represented  by  rQIo  and  the  transmitted 
beam  by  rl.  The  coefficients  r  and  rQ,  accounting  for  the 
cell  losses,  need  not  be  equal. 

The  measured  transmission  of  the  cell-solution  system  is 

T  =  rl/r  I  while  the  true  transmission  is  T  =  e  E  = 

m  oo 

I/IQ  .  Equation  (71)  now  becomes 

JlnT  =  Jin  (-§— )  ~  e(v)cb  (80) 

m  ro 
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If  one  makes  a  series  of  measurements  of  T  for  various  cell 

m 

thicknesses  b,  holding  the  concentration  c  constant  and 
plots  in  Tm  vs  b,  the  resulting  curve  is  a  straight  line 
whose  slope  is  -e(v)c.  The  molar  extinction  coefficient 
e(v)  can  then  be  extracted  from  the  data  by,  viz.  a  least 
squares  fit.  The  ordinate  intercept  of  this  curve  gives 
r/rQ,  but  is  of  little  practical  significance.  As  indicated 
above,  we  have  used  this  method  to  obtain  e(v)  for  tri-n- 
butyl  phosphate.  As  a  check  on  the  validity  of  the  method, 
we  determined  the  molar  extinction  coefficient  e (v)  and  from 
it,  the  extinction  coefficient  k(v)  =  (c/4tt)  e(v)/v  [see 
Eqn.  (72)]  for  pure,  spectral  quality  chloroform  (CHCl^) . 

Our  results  at  3020  cm- 1 (3 . 31y ) [k (3020)  =  0.011]  agree  very 
well  with  those  reported  by  Hauranek  and  Jones  [Spectrochim. 
Acta  32A,  111  (1976)]  [k(3020)  =  0.0115]  using  a  considerably 

more  involved  procedure. 

3.  Other  Corrections 

Using  detailed  computer  calculations  for  measured  and 
simulated  data,  Jones  and  coworkers*  have  considered  fac¬ 
tors  contributing  to  spectrophotometric  errors  in  determin¬ 
ing  the  optical  constants  of  materials  confined  to  thin 
cells  (l-100u).  For  information  only,  we  list  these  factors 
with  our  parenthetical  comments  where  appropriate. 

(a)  Construction  of  precision  (thin)  cells. 

(b)  Distortion  caused  by  interference  and  reflection 
within  the  cell  and  at  cell  surfaces. 

(c)  Polarization  discrimination  of  monochromator. 

(d)  Distortion  caused  by  wedge-shaped  cell,  i.e.,  non¬ 
parallel  windows. 


* 

See  Spectrochima  Acta  32A,  75-123  (1976) 
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(e)  Non-normal  beam  incidence  due  to  beam  convergence  - 
(effects  interference  fringe  pattern  when  b>100y) . 

(f)  Surface  roughness  of  window  faces  -  negligible  in 
mid  and  far-infrared,  i.e.  for  A  -  2.5y  (cell  win¬ 
dows  should  be  demounted  and  repolished  -  every 
three  months ) . 

(g)  Indeterminacy  in  cell  thickness  value  caused  by 
cold  flowing  of  window  material  and  (e)  (cell 
spacing  should  be  measured  before  and  after 
spectrum  is  obtained) . 

(h)  Dispersion  distortion  of  absorbing  solutes  in  non¬ 
absorbing  solvents  [important  for  high  concentra¬ 
tions  (i.e.  c  ~  50  mol%)  in  thin  cells  (i.e. 

b  -  ly)]. 

(i)  Various  computational  errors  generated,  for  ex¬ 
ample,  by  truncation  of  wave  number  interval  in 
computer  program  (not  a  factor  in  this  work) . 

Many  of  the  above-listed  factors  are  of  significant 
importance  only  if  extremely  accurate  (i.e.  archival)  de¬ 
terminations  of  optical  constants  are  sought.  In  these 
cases  considerable  experimental  control  and  computer  time 
are  required. 

B.  Solid  Films 

1.  Reflection  Losses 

Since  free  standing  solid  films  have  only  two  inter¬ 
faces,  as  compared  to  four  in  a  liquid  film  cell,  the  opti¬ 
cal  problems  involved  are  somewhat  simpler.  Reflections 
will  still  occur  at  the  first  and  second  surfaces  and  inter¬ 
ference  fringes  may  appear  if  the  film  thickness  has  the 
appropriate  value.  Fringes  may  be  eliminated  by  (a)  casting 
a  film  directly  on  a  transparent  window  or  (b)  coating  the 
film  with  Nujol  (mineral  oil)  and  pressing  it  against  a 
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thick  window  [Latinski,  C.  ,  Anal.  Chem.  3(),  2071  (1958)  J. 

Of  course  if  the  film  is  thick  enough,  the  fringe  separa¬ 
tion  is  less  than  the  monochromator  resolution,  and  inter¬ 
ference  fringes  are  not  observed. 

Reflection  losses  at  the  air/film  interfaces  may  be 
treated  in  a  manner  similar  to  that  described  above  for 
liquid  film  measurements.  That  is,  the  spectrum  is  scanned 
repeatedly  using  films  of  varying  thickness.  One  then  uses 
the  equation 

’  « 

Jin  Tm  =  £n  -  a  (v)  b  (81) 

to  extract  the  absorption  coefficient  a(v).  The  extinction 
coefficient  k(v)  =  Im[n*(v)]  may  then  be  obtained  [see  Eqn. 

(70) ] ,  if  so  desired. 

2 .  Other  Errors 

The  discussion  presented  above  in  paragraphs  A. 2  and 
A. 3  are  generally  applicable  for  thin,  free-standing  films, 
with  some  minor  modifications  possibly  required.  Therefore, 
they  will  not  be  further  addressed  here.  Of  greater  impor¬ 
tance  here  are  two  experimental  difficulties  related  to  sample 
preparation.  First:  Is  the  (solid)  sample  available  in  thin 
film  form?  Here  "thin"  means  sufficiently  transparent  in  the 
wavelength  region  of  interest  so  that  photometric  errors 
(see  paragraph  A.l)  are  minimized,  i.e.  so  that  T  *  30-40%. 
Second:  Is  the  thin  film  of  uniform  thickness? 

If  the  solid  material  is  insoluble,  such  as  with  some  poly¬ 
mers,  films  may  be  produced  by  rolling  or  pressing  the  bulk 
sample  in  a  molten  state.  Both  methods  offer  positive  and 
negative  features.  For  example,  hot  pressing  of  polymer 
beads  on  powder  may  not  offer  sufficient  temperature  control, 
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allowing  for  thermal  degradation,  or  may  allow  air  pockets 
or  bubbles  to  form  within  the  film,  causing  photometric 
inaccuracies.  Rolling  samples  may  avoid  these  problems, 
but  introduces  orientation  effects,  possibly  undesirable. 
Sectioning  thin  slices  is  possible  if  a  microtome  is  avail¬ 
able  and  if  the  material  is  rigid  enough. 

For  soluble  materials,  casting  films  from  solutions  is  a 
practical  method.  Numerous  techniques  have  been  developed 
to  produce  uniformly  thick  films,  many  of  which  have  been 
reviewed  by  Mano  and  Durao  [J.  Chem.  Ed.,  5j),  228  (1973)]. 

The  choice  of  solvent,  environment  of  cast  solution  and 
condition  of  the  substrate  surface  are  among  the  relevant 
factors  to  be  considered  here. 

III.  EXPERIMENTAL  RESULTS 

A.  Polymethylmethacrylate  (PMMA) 

PMMA  samples  were  obtained  from  cast  sheets  of  Lucite 
supplied  by  the  Plastic  Products  and  Resins  Department  of 
E.  I.  DuPont  de  Nemours  &  Co.  Lucite  pieces  were  dissolved 
in  reagent  grade  toluene.  The  resultant  solution  was  poured 
into  teflon  casting  dishes  which  were  set  upon  a  leveled 
microflat.  The  rate  of  solvent  evaporation  was  controlled 
by  enclosing  the  casting  dishes  within  a  large  glass  con¬ 
tainer  (inverted  fish  tank) .  Sample  thickness  was  controlled 
by  varying  the  volume  of  solution  into  the  casting  dish. 

In  spite  of  the  above  precautions  undertaken  to  obtain 
uniform  films,  the  resulting  sample  thicknesses  exhibited 
considerable  variation.  The  three  PMMA  films  used  to  deter¬ 
mine  a(v)  were  measured  and  found  to  have  the  following 
values . 


b 


0.375  ±  .029  mil  =  9.53 
0.225  ±  .050  mil  =  5.72 
0.169  ±  .067  mil  =  4.29 


0 .74pm  (  8%) 
1.27pm  (22%) 
1.70pm  (40%) 


Film  thickness  was  measured  with  a  Starrett  Gauge  (No.  25- 
209) .  The  average  of  b  "seen  by"  the  optical  beam  is  listed 
above  together  with  the  (%)  deviation.  This  magnitude  of 
the  deviation  in  film  thickness  is  the  principal  source  of 
error  in  the  determination  of  a(v)  for  PMMA  and  is  a  conse¬ 
quence  of  the  need  to  have  films  thin  enough  so  that  suffi¬ 
cient  radiant  energy  is  transmitted  and  reaches  the  detector 
of  the  spectrophotometer  (see  paragraph  II) .  Since  strongly 
absorbing  materials  are  of  primary  interest  in  this  project, 
the  problems  associated  with  producing  "uniform  and  thin" 
films  cannot  be  avoided. 

The  infrared  spectra  of  the  three  films  of  PMMA  were  re¬ 
corded  from  4000  to  800  cm  ^  (2 . 5  -  12 . 5pm)  .  The  value  of 
a(v)  was  calculated  from  Eqn. (81)  at  20  cm  ^  intervals  from 
4000  to  2000  cm  *  (2.5  -5pm)  and  at  10  cm  1  intervals  from 
2000  to  850  cm  ^  ( 5 . 50  -  12 . 5pm)  All  values  reported  were  cal¬ 
culated  by  a  linear  regression  analysis  on  the  data  points 
(£n  T  vs  b) .  The  results  are  presented  in  Fig.  4a  and  4b. 

B.  Tri-n-Butyl  Phosphate 

Solutions  of  tri-n-butyl  phosphate  (TBP)  were  prepared 
by  mixing  purified-grade  solute  with  spectral  quality  sol¬ 
vents  (CCl,^  and  CS2).  Various  concentrations  were  produced 
by  measuring  known  quantities  of  TBP,  using  Van  Slyke  pipets , 
into  calibrated  volumetric  flasks  and  filling  to  level  with 
solvent.  Calibrated,  sealed  liquid  cells  were  used  for  path 
lengths  down  to  approximately  0.1mm.  At  higher  concentrations, 
demountable  cells,  using  various  thinner  spacers,  were  re¬ 
quired.  Where  possible,  path  lengths  were  determined  by 
measuring  cell  fringe  spacing.  For  very  thin  cells,  for 
which  the  fringe  method  is  not  applicable,  the  spacer  thick- 
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ness  was  determined  as  follows.  For  a  given  concentration, 
a  weak  TBP  absorption  band,  adjacent  to  the  band  of  interest, 
is  used  as  an  internal  calibration  standard.  That  is,  its 
absorption  vs.  path  length  behavior  is  known  from  previous 
measurements  utilizing  fixed,  sealed  cells.  The  measured 
absorbance  of  this  weak  band  is  then  used  to  determine  the 
path  length  of  the  demountable  cell.  This  method  has  an  ad¬ 
vantage  in  that  both  thickness  of  cell  and  absorbance  of  de¬ 
sired  band  are  measured  simultaneously  and  in  a  relatively 
short  time  interval.  Thus,  errors  introduced  by  attempting 
to  duplicate  a  pressure  required  to  produce  a  given  cell 
spacing,  and  cell  distortion  are  minimized.  We  have  found 
this  method  to  be  quite  satisfactory  for  very  thin  demount¬ 
able  cells  and  have  used  it  for  spacings  to  .0013mm.  As  re¬ 
commended  by  Hawranek1  no  compensating  solvent-containing 
cell  was  used  in  the  reference  beam. 

Figure  5(a)  shows  the  infrared  spectrum  of  a  capillary 
film  of  pure  TBP.  The  absorption  bands  are,  as  expected, 
characteristic  of  the  butyl  [(C4H9)  ]  and  phosphate  [(P04)-] 

groups.  Consistent  with  our  interest  in  strongly  absorbing 
compounds,  we  have  selected  the  three  most  intense  bands  in 
the  2-5  and  8  -  12pm  regions  for  detailed  analysis.  These 
are  at  2950cm  1280cm  ^  and  1030cm  ^  and  are  assigned  to 

the  CH,  P  =  0  and  P  -  0  groups  respectively. 

In  studying  the  absorptivity  of  a  possible  obscuring 
agent,  concentration  effects  are  of  major  interest,  both  from 
a  scientific  and  practical  point  of  view.  Therefore,  we  have 
selected  a  range  of  concentrations  from  very  dilute  (.0073  M) 
to  pure  liquid  (3.65  M)  for  investigation. 

For  each  solution  of  TBP,  infrared  spectra  were  recorded 
using  at  least  three  liquid  cells  with  different  path  lengths. 
Multiple  scans  were  obtained,  and  average  measured  transmit¬ 
tance  values  for  each  of  the  three  bands  of  interest,  were 
calculated.  Eqn. (80)  was  used  to  obtain  e (v)  for  each  con- 
1 [Spectrochimia  Acta  32A,  75  (1976) ] 
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Fig.  5a  -  Infrared  Spectrum  of  a  Capillary  Film  of  Tri-n-Butyl 
Phosphate.  (Bands  used  for  analysis  are  indicated.) 


centration,  using  a  least  squares  fit  of  the  Jin  TM(v)  vs.  b 
data.  The  results  are  shown  in  Figure  5(b).  Systematic 
errors  were  estimated  from  the  uncertainties  in  transmit¬ 
tance,  path  length  and  concentration.  Of  these,  the  uncer¬ 
tainty  in  path  length  b  was  the  dominant  contribution.  The 
fractional  determinate  error  in  e(v)  was  approximately  4-6% 
for  the  bands  studied. 

As  discussed  above  in  paragraph  II,  in  the  Bouguer-Beer 
law,  Eqn.(71),  the  value  of  the  molar  extinction  coefficient 
should,  ideally  be  independent  of  concentration.  In  real 
systems,  however,  one  expects  deviations  from  ideal  behavior. 
Figure  5(b)  suggests  a  behavior  where  the  molar  extinction 
coefficient  is  constant  at  low  concentrations  (normal  Beer’s 
law  behavior) ,  undergoes  a  transition  to  lower  values  as  the 
concentration  increases  and  "saturates"  at  high  concentrations. 
Note  that  the  transition  regions  for  the  1030cm  and  1280cm  , 

obtained  using  the  same  solvent,  are  very  similar  in  shape. 

The  onset  of  the  transition  region  for  the  2950cm  ^  band, 
however,  occurs  at  a  considerably  higher  concentration.  Note 
also  that  for  this  band  the  solvent  used  was  CC1  while  CS, 
was  used  for  the  two  lower  frequency  bands.  Both  CC14  and  CS2 
are  considered  to  be  non-polar  solvents.  If  dipolar  inter¬ 
actions  between  solute  and  solvent  are  eliminated,  we  must 
then  consider  other  factors  responsible  for  the  observed 
behavior. 

In  measurements  of  the  OH  stretching  band  in  n-butanol 
dissolved  in  various  non-polar  solvents,  Graja  and  Malecki 
[Chem.  Phys.  Lett.  _3i'  373  (1975)  ]  found  the  molar  extinc¬ 
tion  coefficient  versus  concentration  curve  to  be  charac¬ 
terized  by  two  distinct  regions.  At  low  concentrations 
£ (v0H)  was  independent  of  concentration.  At  higher  concen¬ 
trations,  e  (v^)  decreased  smoothly  with  concentration. 

These  investigators  attributed  the  behavior  to  self-associa¬ 
tion  of  the  solute  molecules.  In  effect,  at  low  concentra- 
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tions  the  OH  group  is  quasi-free  so  that  its  dipole  moment 
oscillation  is  unhindered.  At  high  solute  concentrations, 
tie  n-Bu  OH  molecules  interact  and  restrict  the  OH  stret¬ 
ching  motion  which  in  turn  causes  a  reduction  in  the  molar 
extinction  coefficient. 

In  Figure  5(b),  we  can  clearly  see  the  presence  of  two 
distinct  regions,  at  high  and  low  concentrations,  separated 
by  a  transition  region.  For  example,  the  2950cm  band  has 
a  molar  extinction  coefficient  of  approximately  1365  5,/mole-cm 
in  dilute  solution,  and  a  value  of  ~  950  £/mole-cm  as  a  pure 
liquid.  Thus,  the  extinction  due  to  the  CH  stretching  modes 
in  TBP  decreases  by  a  factor  of  about  three  over  the  range 
of  concentrations  studied.  The  similar  behavior  of  the 
1030cm  ^  and  1280cm  ^  bands  indicate  that  the  motions  of  the 
phosphate  ion  modes  in  TBP  are  also  dependent  upon  concen¬ 
tration.  The  difference  in  the  curves,  other  than  the  mag¬ 
nitudes  of  e(v)  which  are  determined  by  electric  dipole 
transition  probabilities,  is  the  point  at  which  the  transi¬ 
tion  region  begins.  For  the  1030cm  1  and  1280cm  1  bands, 
the  transition  begins  at  a  concentration  of  about  .03 M while 
for  the  2950cm  ^  band  the  corresponding  concentration  is 
about  0.3  M.  We  believe  this  difference  is  due  to  a  solvent 
effect  related  to  the  shapes  of  the  respective  solvent  mole¬ 
cules.  Carbon  tetrachloride,  a  tetrahedrally  bonded  struc¬ 
ture,  can  roughly  be  considered  as  a  sphere  of  diameter 

O 

~3.6A.  Carbon  disulphide,  also  non-polar,  is  a  linear  mole- 

O 

cule  with  an  axial  length  of  about  3.6A.  It  would  seem  that 
the  packing  of  CS?  "rods"  between  TBP  molecules  can  be  ac¬ 
complished  such  that  solute  molecules  can  interact  at  a 
lower  concentration  than  possible  for  the  CC1  spheres. 

4 

Viewed  alternatively,  the  CC14  molecules  are  more  effective 
at  preventing  TBP  'association'  and  consequently  the  transi¬ 
tion  to  the  lower  extinction  state  is  held  off  to  higher  con¬ 
centrations  . 
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The  solvent  effects  described  here  are  manifest  in  the 
experimental  procedure  to  determine  the  material  parameter, 
e (v) .  As  such  they  have  limited  applicability  to  possible 
field  use  of  TBP  as  an  obscurant.  However,  the  fact  that 
the  extinction  coefficient  is  concentration-dependent  is  of 
interest.  In  particular,  these  results  show  that  it  can  be 
erroneous  to  evaluate  the  molar  extinction  coefficient  by 
extrapolating  the  curves  of  e(v)  versus  concentration  to 
dilute  conditions.  Also,  in  the  case  of  TBP,  'more  is  not 
better ' . 

C.  Cellulose 

Samples  of  regenerated  cellulose  film  were  obtained  from 
the  Film  and  Packaging  Division  of  FMC  Corporation.  Two  film 
thicknesses  were  specially  prepared  for  our  use  and  contained 
no  plasticizer.  The  films  were  dried  under  tension  in  FMC 
Corporation’s  Laboratory.  While  the  moisture  content  of  cel¬ 
lulose  is  very  sensitive  to  atmosphere  humidity,  we  were  un¬ 
able  to  determine  any  appreciable  difference  between  the 
spectrum  of  a  film  stored  in  air  and  one  which  had  been 
dryed  for  24  hours  under  vacuum.  No  additional  drying  at 
elevated  temperatures  was  attempted  by  us. 

Cellulose  is  a  naturally  occurring  material  having  a 
typical  fibre  structure.  It  is  the  main  component  of  the 
cell  walls  of  most  plants  and,  thus,  has  been  estimated  to 
comprise  about  33%  of  all  vegetable  matter.  Cellulose,  to¬ 
gether  with  starch  and  glycogen  are  polysaccharides  composed 
of  glucose  only.  The  glucose  units  are  linked  by  6-glyco¬ 
side  bonds  resulting  in  a  repeat  or  monomer  unit  represented 
by 
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and  referred  to  as  a  cellobiose  unit. 

Industrial  cellulose,  as  used  here,  is  extracted  from 
wood  pulp  by  a  straight-forward,  but  carefully  controlled  pro¬ 
cess.  Our  laboratory  is  not  equipped  to  perform  the  extrac¬ 
tion  procedure.  The  films  kindly  prepared  for  us  by  FMC 
Corporation  were  received  with  average  thicknesses  of  0.45 
mil.  (.0011  cm)  and  0.95  mil.  (.0024  cm).  Normally,  such 
optical  path  lengths  for  transmission  measurements  are  suit¬ 
able  for  determining  absorption  coefficients.  However,  cel¬ 
lulose  is  highly  absorbing  in  the  3600-3000  cm-1  (2.8-3.33y) 
and  1200-950  cm-1  (8.3-10.5y)  regions.  This  is,  of  course, 
a  highly  desirable  property  in  terms  of  the  goals  of  the 
overall  project.  The  low  transmission  of  the  samples  in  these 
regions  made  quantitative  measurements  impractical.  Therefore, 
in  Figs.  6  a  and  6  b,  where  the  behavior  of  ot(v)  for 
cellulose  is  presented,  those  regions  for  which  numerical 
values  are  not  reported  are  indicated  by  dashed  lines.  It 
may  be  possible  to  produce  samples  of  cellulose  that  are  thin 
enough  to  yield  transmissions  of  30%  in  those  highly  absorb¬ 
ing  regions,  yet  are  mechanically  stable  enough  to  be  treated 
as  free  films.  However,  to  do  so  should  require  a  special 
order  to  FMC  Corporation  without  additionally  imposing  on 
their  generosity  and  cooperation. 

D.  Polyvinyl  Alcohol 

Polyvinyl  alcohol  (PVA)  resin,  in  powdered  form,  was  do¬ 
nated  by  American  Hoechst  Corp.  Material  received  was  de¬ 
signated  as  Moviol  4-88  and  characterized  by  a  viscosity  of 
4  cP  (centipoise) ,  for  a  4%  aqueous  solution  at  20°C,  and  an 
88%  degree  of  hydrolysis.  PVA  solutions  were  produced  by 
adding  the  powder  to  stirred  cold  water,  then  heating  the 
mixture  to  85-90°C  until  dissolution  was  complete.  The  re¬ 
sulting  solution  is  then  poured  into  casting  dishes  upon 
which  films  are  formed  when  the  solvent  has  evaporated. 
Approximately  1%  solutions  (by  weight)  were  used.  Sample 
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thickness  was  controlled  by  varying  the  volume  of  solution 
poured  into  the  casting  dish. 

Three  films  were  produced  with  thicknesses  of  23.6±2.5ym, 
3.81±0.00iJm,  and  2. 90  +  0. 61ym.  The  errors  in  the  thicknesses 
were  determined  by  measuring  only  the  variations  of  those 
areas  of  the  films  "seen"  by  the  optical  beam  of  the  in¬ 
strument.  The  optical  density  of  the  23.6um  film  was  too 
large  to  allow  for  quantitative  results  to  be  reported  over 
the  wavelength  ranges  of  interest.  Therefore,  values  of 
a(v)  reported  in  Fig.  7  a  and  b  represent  data  obtained  from 
the  two  thinner  films. 

E.  Ureaformaldehyde 

Powdered  ureaformaldehyde  (UF)  resin  was  supplied  by 
R.  W.  Doherty  (CSL-Edgewood  Arsenal)  and  characterized  as 
having  an  approximate  particle  size  of  4ym  diameter.  The 
resin  is  crosslinked*  and  so  is  not  soluble.  Because  of  its 
relatively  poor  characterization,  ureaformaldehyde  could  only 
be  investigated  as  a  powder  suspended  in  a  matrix  (Nujol  Oil) . 
Quantitative  measurement  techniques  for  such  powders  are 
presently  being  developed  in  our  laboratory.  Consequently, 
only  a  qualitative  spectrum  of  UF  is  presented  (see  Fig. (8a  and  8b) 
at  this  time. 


^Personal  communication  with  J.  Schwendeman,  Monsanto  Corp. 
Dayton ,  Ohio 
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Fig.  7  a  -  Absorption  Coefficient  of 
Polyvinyl  alcohol.  2.8-4.5ym 
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Fig.  7b  -  Absorption  Coeffici¬ 
ent  of  Polyvinyl  alcohol. 
6-13um _ _ 
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Fig.  8b  -  Infrared  Spectrum  of 
Ureaformaldehyde .  Sample 
prepared  as  a  KBr  pellet. 


GLOSSARY  OF  SELECTED  TERMS 


Symbol 

Quantity 

Units 

Definition 

X 

Wavelength 

pm  =  micrometers 

V 

Wavenumber 

cm-1 

v  = 

1 

X 

c 

Speed  of  light 

meters 

sec 

C  = 

1//  e  p  = 

0  0 

3x10s  meters/s 

f 

Frequency 

Hertz,  Hz 

f  = 

C 

~ir 

Id 

Angular  Frequency 

radians 

sec 

(d  = 

2vf  =  2ttCv 

K 

Wave  Vector, 
Magnitude 

radians 

cm 

K  = 

2tt 

X 

K*  (w) 

Wave  Vector, 
(Complex) 

radians 

cm 

K*  = 

b>  * 

=  —  n* 

n 

Index  of 

Refraction 

none 

»4 

V  =  speed  of 
light  in  medium 

n*  (v) 

Index  of  Refrac¬ 
tion  (Complex) 

none 

n*(v)  =  n(v)+ik(v) 

k(v) 

Extinction 

Coefficient 

none 

k  (v) 

=  Im[n* (v)  ] 

A 

Absorbance  or 
Optical  Density 

none 

A  = 

1Og10  ~ T 

T 

Transmission 

none 

T  = 

I  / ratio  of  \ 

I  /  transmitted! 
to  incident! 

\ power  J 

ot  (v) 

Absorption 

Coefficient 

cm"^ 

a  (v) 

=|[a(v)=ijk(v) 

'  \ 

=  4irvk  (v)  j 
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Symbol 

Quantity 

Units 

Definition 

b 

Optical  Path 
Length 

cm 

e  (v) 

Molar 

Extinction 

Coefficient 

liter 

mole-cm 

,  »  A 

E(v)  =  be 

fe  (v)  •  ^Lyk  (v 

c 

Concentration 

mole 

liter 

1 

A 


Appendix  A:  Development  of  General  Wave  Equation  to  Describe 
The  Propagation  of  Light  in  Dielectric  or 
Conductive  Media. 


The  behavior  of  electromagnetic  waves  in  any  medium  is  des¬ 
cribed  by  the  classical  wave  equation,  which  in  turn  is  derived 
from  Maxwell's  Equations.  In  this  section  the  development  of 
the  appropriate  wave  equation,  of  interest  to  this  work,  will 
be  briefly  presented. 

The  material  parameters  used  to  characterize  the  medium 
with  respect  to  electromagnetic  phenomena  are,  in  SI  units,  the 
relative  permitivity  =  e/eo(eo=  8.85  x  10  12  farad  m-1  = 

permitivity  of  free  space)  and  the  relative  permeability 
=  u/uo(uo=  1.26  x  10-6  henry  m” 1  =  permeability  of  free 
space).  The  velocity  of  propagation  of  the  wave  in  this  medium 
is 

V  =  --1 -  =  1  (A- 1 ) 

/ey  Vkzk  y 
e  o  m  o 

which,  for  free  space  (in  vacuo  K  =  k  =1)  reduces  to 

e  m 


-  3  x  10 8  m/s 


The  index  of  refraction  of  the  medium  is  given  by 


(A-2 ) 


u  = 


C 

V 


K 

e  m 


For  most  transparent  optical  media  K^=  1  so  that 


n  -  -s/e/co 


(A-3 ) 


(A-4) 


This  result  is  generally  only  valid  for  dilute  gases  and  non¬ 
polar  solids,  where  high  static  polarizability  is  present.  In 
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addition,  n  is  a  function  of  frequency  (dispersion)  so  that  to 
characterize  the  medium  optically,  these  effects  must  be  account¬ 
ed  for  and  (A-4)  is  seen  to  be  a  special  limiting  case. 

Maxwell's  Equations  for  a  non-magnetic ,  electrically  neu¬ 
tral  medium  can  be  reduced  to 

V  x  £  =  -  UQ  -|-S-  (A-5a) 

7  x  »  =  eo  TT"  +  TT"  +  3  (A_5b) 

V  •  E  - - 1 -  7  •  P  (A-5c) 

eo 

7  •  H  =  0  (A-5d) 

where  E  is  the  electric  field  intensity,  H  is  the  magnetic  field 
intensity, 

£  =  ( e— e Q )  E  =  xeQE  (A-6) 


is  the  electric  polarization  (x  =  — | —  -  1  =  electric  suscepta- 
bility)  and 

3  =  oE  (A- 7) 

is  the  current  density  (a  =  electrical  conductivity) . 

By  manipulating  the  above  equations  [take  the  curl  of  (A-5a) 
and  time  derivative  of  (A- 5b)  and  eliminate  H  between  the  two 
resulting  equations]  one  can  get 


7  x 


(7xE)  +  — 
C2 


a  2e 

3  t2 


(A- 8) 


From  the  vector  identity 


V  x  (V  x  E)  =  7  (7 *E)  -  72  E 


(A-9) 
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and  the  fact  that  V*E  =  0  for  an  electrically  neutral  media,  it 
follows  that 


V2E 


1  3  2E 

C2  3  t2 


3  2P 
3  t2 


+  U. 


± 1 
3  t 


(A-10) 


The  two  terms  on  the  right  hand  side  of  (A-10)  are  called 
source  terms  and  are  interpreted  as  follows: 


(a)  For  a  non-conducting  medium,  -  accounts  for 

^  A.  2 

the  effect  of  polarization  charges  and  is  the 

important  source  term,  being  responsible  for  optical 
effects  such  as  dispersion,  absorption,  etc. 

3  T 

(b)  For  a  conducting  medium,  -y£  is  the  important  term 
and  is  responsible  for  the  high  reflectivity  and  optical 
opacity  of  metals. 

(c)  For  semiconductors,  both  source  terms  must  be  considered. 

Equation  (A-10)  is  the  starting  point  for  discussions  of  dis¬ 
persion  and  metallic  absorption  considered  in  this  work. 
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